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Abstract 

A new class of time-symmetric solutions to the initial value constraints of 
vacuum General Relativity is introduced. These data are globally regular, 
asymptotically flat (with possibly several asymptotic ends) and in general 
have no isometries, but a U(l) x U(l) group of conformal isometries. After 
decomposing the Lichnerowicz conformal factor in a double Fourier series on 
the group orbits, the solutions are given in terms of a countable family of 
uncoupled ODEs on the orbit space. 
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If one casts General Relativity in canonical form one finds that the Einstein equations 
split into two parts. One part, the constraints, has no time derivatives in it: the equations 
are conditions on the initial data. The other part, the evolution equations, determine how 
this initial data develop in time. The constraints are a highly underdetermined system which 
can be turned into elliptic equations by a judicious choice of the free data. 

The simplest version of the constraints occurs in the so-called moment-of-time-symmetry 
(or time-symmetric) problem, where the velocity part of the data, namely the extrinsic 
curvature of the initial 3- manifold M', is set to zero and the only remaining equation is 
R[g] = 0, with g a Riemannian metric on the Cauchy slice and R the scalar curvature of g. 

The standard approach, pioneered by Lichnerowicz, is to select a base metric g and 
construct a solution metric g' by means of a conformal transformation g' = ip^g. One then 
gets a linear elliptic PDE, namely 

L^:=(-A,+i%])^ = 0, ^>0 (1) 

for the conformal factor. In the asymptotically flat case studied here one in addition assumes 
that g' tends to the flat metric at infinity. One then requires that ip, in addition to obeying 
Eq. ([[]), go to one in that limit. 

Simplification should occur if one imposes some symmetry on g. Completely reducing 
this 3 dimensional PDE problem to a one dimensional ODE by imposing a symmetry with 
2 dimensional orbits does not really work. In the case of spherical symmetry this would, by 
the Birkhoff theorem, leave us with data for the Schwarzschild solution. In the cylindrically 
symmetric case, using the boundary condition, only Minkowski space remains [|IJ. Thus, in 
the standard approach, the simplest solution that has been identified is to assume that the 
free data is axially symmetric which reduces Eq. ([[]) to a two dimensional PDE, and much 
effort, especially by the numerical community, has been devoted to solving it. 

In this paper we wish to introduce a new symmetry in the base metric which gives us 
a large class of non-trivial solutions while at the same time allowing us to benefit from 
the underlying linearity of the equation. With this symmetry we can reduce the equation 
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to a set of uncoupled linear ODEs. This is a remarkable simplification which offers major 
benefits in a whole range of problems. The key idea is to use the freedom in choosing the 
base metric g in a conformal equivalence class, particularly one can choose a base metric 
that lives on a compact manifold which can be opened up to an asymptotically flat one by 
conformal decompactification (analogous to stereographic projection). 

The physical data we are considering are data defined on a 3-dimensional manifold M' 
which may have several ends and are asymptotically flat near each of them. They possess two 
commuting, orthogonal conformal Killing vector fields. Their action extends to a U(l) x U(l) 
action on the many-point compactification S 3 of M' . After a suitable conformal rescaling 
of the metric, this action becomes in fact an isometry, and is thus of the "polarized Gowdy 
type" @. In the special case of just one point-at-infinity, which is a fixed point of one 
U{1) factor, these data have this U(l) factor as an isometry and are thus a subclass of Brill 
waves. The fact that the group of all conformal symmetries of the physical data extends to 
an isometry group only after conformal compactification turns out to be a blessing rather 
than a drawback. 

Even in the absence of symmetry it is useful to replace Eq. (J!]), on R 3 say, together with 
the boundary condition if) — > 1 at infinity, by the equation 

L g G = 4vt(5a, G > (2) 

on the one-point compactification S 3 of M' with g now being a Riemannian metric on S 3 and 
5 \ the Dirac delta function concentrated at some point A e S 3 . It is then easy to see, using 
the known asymptotic behaviour of G near A, that the metric g' ab = G A g a b on R 3 = S 3 \A is 
asymptotically flat near A and satisfies R[g'] = on account of Eq. (fj). A solution G > is 
known to exist, iff Ai, the lowest eigenvalue of L g , is positive ||. In particular, this existence 
is insensitive to the choice of the point A which we will use later when we construct solutions 
with many asymptotic ends. 

To solve Eq. (|2|) on the compact manifold S 3 we subtract from G a quantity Gi oc , taken 
to be a suitable approximation to the Hadamard fundamental solution Q of L g near A, 
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arbitrarily extended to all of S 3 . By including a sufficient number of terms in Gi oc , we can 
arrange for the r.h.s. in 

L g (p = L g (G - doc) = Pa (3) 

to be as smooth at A as we please. Solving the regular equation @ on S 3 numerically the 
usual problems arising from having to use a finite grid for an infinite region (R 3 in this case) 
are avoided. Moreover, in the presence of a U(l) X U(l) isometry we can introduce adapted 
coordinates, which will allow us to reduce the PDE (j3|) to a system of uncoupled ODEs. 

For the basic setup it is best to start by obtaining the standard metric g ab on S 3 by inverse 
stereographic projection from R 3 with flat metric g' ab . Let z,r,(p be standard cylindrical 
coordinates on (R 3 ,g'). Taking = 4(1 + 4(r 2 + z 2 ))' 1 , the metric g ab = Q 2 g' ab extends to 
the standard one on S 3 . Furthermore if = (J^) a extends to a Killing vector on (S 3 ,g ab ). 
With p 2 = n a r) b g ab we have that rj a p a = and 

DaVb = Zp^PlaVb], Pab = P^VaVb ~ P9abi ( 4 ) 

where we use the definitions 

O q o o o 

D c g a b = 0) Pa = DaP, Pab = D a DbP- 

Next define 

r = P~ 10 e abc PbVc (5) 

obeying = 1 ~ p 2 =: cr 2 , ?Pa = 0, = 0, and 

Da£ 6 = 2 ( 7-V [o £ 6] . (6) 

It thus follows from equations (TJ,|^^) that (£, rj) form a pair of commuting, orthogonal 
and hypersurf ace-orthogonal Killing vectors of g ab spanning p = p = const, where po ranges 
between and 1. For po 7^ 0, 1 these are flat tori of constant mean curvature changing sign 
at the Clifford torus for which p = 775 • The sets p = (p = 1) are closed lines on which 
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n a (£ a ) is zero (the "axes"). In fact, p = (respectively p = 1) are linked great circles on 
S 3 , corresponding to the z-axis (respectively the circle z — 0,r — 1/2) after stereographic 
projection. Using coordinates (p, x, V 9 ) where £ a = {-§^) a and x = const, are orthogonal to 
£ a we find that Xi like <p, has periodicity 2tt. The coordinate x can be fixed by requiring 
A, the point at infinity under stereographic projection, to be at p = 0, x — 0. The totally 
geodesic surfaces x = const, (<p = const.) are metric hemispheres intersecting at p = 1 
(p = 0) with {x = en} U {x = ol + 7r} ({<p = f3} U {<p = /? + tc}) being smoothly embedded 
S^'s. Under stereographic projection the surfaces x — const, do not change topology except 
for {x = 0} U {x = tt} which gets decompactified into {z = 0, r > 1/2} U {2; = 0, r < 1/2}, 
i.e. the equatorial plane z — 0. 

In the (p, x, <p) coordinates the metric takes the form 

= + I 1 ~ P )XaXfe + P <Pa<P&- (7) 

We now generalize Eq. (^) by setting 

+ p 2 <p <P6 (8) 

where A is a constant and q is a smooth function of p 2 with q\ p= o = 0. Eq. (g) then defines 
a smooth metric on S 3 which, again, has a U(l) X C/ (1) symmetry with the same properties 
as above, with the exception that the hemispheres x = const, and <p = const, are not now 
metric half-spheres and that maximal tori p = p occur whenever 

Ap(l - p 2 )g' + 1 - 2p 2 = 0, where g' = ^. 

ap 

The metrics defined by (|j) are the free data we wish to consider. The first step now 
consists of checking whether \\(g) > 0. This is known to be true for A\ < A < A 2 0, 
where A\ < 0, A 2 > are numbers dependig on q, both finite provided q ^ 0. In the present 
case \i(A) is determined by an ODE: by virtue of non- degeneracy of ground states, the 
corresponding eigenfunction has to share the U(l) x £7(l)-symmetry of the problem and is 
a function merely of p. (In contrast, in Eq. (§), the original symmetry is broken by one's 
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choice of A.) It is known that for A sufficiently close to A\ or A2, the constant term in G, i.e. 
the ADM mass gets arbitrarily large and the physical metric g' develops trapped surfaces 
(see Ref. §). 

In the next step we subtract a suitable Gi oc from G and try to solve Eq. (|3|). We can 
now expand and p\ in a Fourier series on the square torus (x € [0, 2ir],(p G [0, 2ir]). Since 
L g commutes with ^ and ^ this results in an ODE on p E [0, 1] for each Fourier mode, 
supplemented by boundary conditions at p = and p = 1 needed to ensure regularity of (f> 
on the two axes. By uniqueness || an isometry of g ab carries over to g' ab iff A is a fixed point 
of this isometry. Thus, if A is in a general position, g' has no isometries whatsoever. On 
the other hand, taking A to be at p = 0, x = 0, say, we keep ^ as a Killing vector, and the 
resulting data are a special case of Brill waves j7|. The vector field gp survives on R 3 = S 3 \A 
merely as a conformal vector field which is incomplete along the axis p = 0. There is however 
the discrete "mirror" isometry (p, Xi f) l— * (p? 27r — x,f), leaving the equator invariant. 

We can generalize the above to the case where there are several asymptotic ends. We 
simply choose a finite number of points Aj (i = 1, . . . , N), solve numerically for the corre- 
sponding Green function Gi and write 




where the "source strengths" q are arbitrary positive numbers. When A = 0, these are 
simply the many-black-hole solutions discussed by Brill and Lindquist ||. 

If all Aj lie on the two axes, each Gi in G = c iG% corresponds to an axially symmetric 
(Brill-type) solution. It is perhaps useful to observe that Brill's simple positive-mass proof 
immediately gives positivity of mass also in this non-axially symmetric, topologically 
nontrivial situation. Another useful observation, in the case of a general distribution of 
Aj's, is that not all of the G^s have to be computed separately: for any A 1; A 2 lying on 
the same orbit of the isometry group, G2 is given in terms of G\ by a product of suitable 
rotations in ip and x- Take e.g. N = 2 with A x at (p = 0, x — 0) and A 2 at (p = 0, x — 
In the special case a = n, the physical metric g' ab has the mirror symmetry x l— > 27r — x- 
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When ci = C2, it has the inversion symmetry {%} U + 7r} i— > {a — x} U { a ~ X + ^ji 
leaving the "throat" {x = ol/2} U{a/2 + 7r} invariant. Thus the "throat" is totally geodesic 
with respect to g' ab , in particular a minimal embedded 2-sphere (a "horizon"). We remark 
that, when A = 0, the above discrete symmetries are present irrespective of C\, c 2 and a. 
Because, then, by combining a homothety with a proper conformal motion coming from the 
Killing vector on (R 3 ,g' ab ), we can always arrange for c[ = c' 2 , a' = ir. The corresponding 
physical data is of course nothing but a time-symmetric slice of the Kruskal spacetime with 
mass m(ci, C2, a) = 2V2c lC2 (l — cos a) ' . When A 7^ our solutions for c\ 7^ C2 are not 
inversion symmetric. Placing sources of equal strength at p = 0, aij = ^p, i = 0, 1, . . . , N—l, 
we find in an analogous manner that on R 3 \ 1J" =1 Aj each asymptotic end is surrounded by 
a throat. Thus, viewing Ai as "infinity" and the other Aj's as "black holes" or "particles" 
we can in particular say that these objects, viewed from infinity, are so close together that, 
in addition to a horizon for each of them individually, there is one surrounding them all. 

Finally we want to place the solutions found here in a more general context. From the 
conformal-compactification viewpoint adopted here, looking for time-symmetric initial data 
with symmetries naturally leads to looking for metrics g on a compact 3-manifold M, not 
necessarily S 3 , with a group Q of conformal symmetries. Suppose (M,g) is not conformally 
diffeomorphic to (S 3 ,g) (in which case the one-point conformal decompactification would 
be data for Minkowski spacetime. Then, by virtue of a theorem due to Obata |J, Q acts by 
isometries on a metric g\ in the conformal class of g. Thus, calling g\ again g, we are left with 
isometries of (M, g). When Q has 3-dimensional orbits, it has to be one of the Bianchi types 
existing on a compact manifold. Whether they too can be usefully employed as background 
metrics for conformal decompactification, we have not investigated. The next case is the 
one of 2-dimensional orbits which has been fully classified [[HJ. Taking M to be S 3 , which 
is the topologically trivial case in the present context, Q can be 5*0(3) or U(l) x U{\). The 
first case is again conformally, whence physically, trivial. The second case, provided the 
two Killing vectors can be chosen to be orthogonal, is the polarized Gowdy case, which is 
conformally diffeomorphic to the one used here. 
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If M is chosen to be S 2 x S 1 and Q = SO (3), (M, g) is conformally equivalent to S 2 x 5(a), 
the standard S 2 times the circle of radius a. This, in physical terms, corresponds to the 



Misner wormhole ||11|| . The case Q — V (1) xU r (l) is a generalization of the Misner wormhole, 
and is similar to the case studied here. 

Finally, when M = T 3 , it follows from a fundamental theorem of Schoen and Yau [T2| 
that Ai(g) < for all metrics g, so that no positive Green function exists. 

Let us summarize: We have presented here a class of initial data for General Relativity 
which in general have no symmetries whatsoever and for which the constraints are equivalent 
to an infinite system of uncoupled ODEs. In forthcoming work by one of us (S.H.), the 
procedure outlined here is carried out for some classes of initial data having the required 
U(l) x U(l) conformal symmetry. This will involve a numerical study of the existence 
question, i.e. the sign of Xi(g), and the actual construction of initial data. Also included will 
be the location of apparent horizons and tests for the isoperimetric inequality of black holes 
and the hoop conjecture. Since our method is fundamentally different from the one usually 
employed in Numerical Relativity based on solving PDEs directly on 2 or 3 dimensional 
grids, we expect our solutions to be useful testbeds for ones obtained by these standard 
methods. Clearly our method so far makes decisive use of the linearity of the Lichnerowicz 
equation (§) in the time-symmetric case. In the future we plan to study the nonlinear case of 
maximal non-time-symmetric initial data having a U(l) x U(l) symmetry of the background 
metric and background extrinsic curvature. Whether the lines of thought followed in the 
present work offer simplifications also for the evolution problem is at present unclear. 

This work was supported by Fonds zur Forderung der wissenschaftlichen Forschung, 
Project No. P09376-PHY. 
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